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Abstract 



In this paper we investigate the survival probability, in high-dimensional statistical physical 
models, where 9n denotes the probability that the model survives up to time n. We prove that if the 
r-point functions scale to those of the canonical measure of super-Brownian motion, and if a certain 
t:;;^- , self-repellence condition is satisfied, then n9n — >■ 2/(AF), where A is the asymptotic expected number 

of particles alive at time n, and V is the vertex factor of the model. Our results apply to spread- 
out lattice trees above 8 dimensions, spread-out oriented percolation above 4 + 1 dimensions, and the 
spread-out contact process above 4-1-1 dimensions. In the case of oriented percolation, this reproves 
a result by the first author, den Hollander and Slade (that was proved using heavy lace expansion 
i^H ' arguments), at the cost of losing explicit error estimates. We further derive several consequences of 

, our result involving the scaling limit of the number of particles alive at time proportional to n. Our 

proofs are based on simple weak convergence arguments. 



1 Introduction and results 

> 

Q"^ I A celebrated result by Kolmogorov [32] states that the probability 9n that a Galton- Watson branching 
■ process with offspring distribution having mean 1 and variance 7, starting from a single initial particle, 
survives until time n satisfies n9n — )-2/7asn— )>oo (see also [39, Theorem H.I.I.]). A related classical 
I result by Yaglom [43] states that the population size Nn at time n is such that, conditionally on survival 
up to time n, the random variable n~^Nn converges weakly to a random variable Y having an exponential 
distribution with mean 7/2. Thus, the probability of survival up to time n decays like 1/n, while on the 
event of survival, the number of alive particles grows proportional to n. In this paper, we study extensions 
of this result, and their ramifications, to general spatial statistical mechanical models in sufficiently high 
dimensions. 

, We next define the scaling limit of the particle numbers for critical Galton- Watson trees. The proba- 

bility of a particle surviving is rather small, and in the literature, two constructions have been investigated 
to resolve this problem. The first construction to deal with the vanishing survival probability is to start 
with a large number of particles, i.e., take A^o = Inx] , where x > 0. In this case, at any time t > 0, the 
number of particles at time whose lineage survives until time t has an approximate Poisson distribution 
with parameter 2x/'y. Then, the process {Ntn/n)t>o converges in distribution to Feller's branching diffu- 
sion [13], which is the unique solution to a stochastic differential equation describing a continuous-state 
branching process (see also [35] for related results). The second construction to deal with the vanishing 
survival probability is to multiply the measure by a factor of n, making sure that the measure of the 
event of survival to time proportional to n converges to a finite and positive limit. Then, the process 
{Ntn/n')t>o converges in distribution, where the notion of convergence in distribution is defined in terms 
of convergence of integrals of bounded continuous functions having support on paths that survive up to 
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time e > 0. The resulting measure is a cr-finite measure rather than a probabiUty measure, and is called 
the canonical measure of the branching process in reference to canonical measures appearing in infinitely 
divisible processes (see e.g. [31]). We can retrieve a probability measure by 'conditioning' the measure on 
surviving up to time 1. 

While the two constructions are quite different, they are closely related. Indeed, in the first construc- 
tion (conditionally upon survival to time 1) take any of the Poisson 2x/7 initial particles whose lineage 
survives until time 1. Then the distribution of its rescaled numbers of descendants is identical to that in 
the canonical measure conditioned to survive up to time 1. 

The models we consider will be spatial. Embedding the branching process into Z'^, with the initial 
particle located at the origin, G Z*^, and where the offspring of any given particle are independently 
located at neighbors of that particle in Z*^, we obtain a branching random walk. Since multiple occupancy 
can occur, the state of this process at time n is best described by a (random) measure, where the measure 
of any subset of is the number of particles of generation n located in that set. With appropriate 
rescaling of space, time, mass (associated to each particle), and of the underlying law, we obtain a 
sequence of finite (no longer probability) measures Watanabe [42] shows that the measures fin 
converge weakly to a measure No on the space of measure valued paths {Xt)t>o that survive for positive 
time, i.e. S = mf{t > 0: Xt{l) = 0} > (where Xt{f) = J fdXt). The measure No is called the canonical 
measure of super- Brownian motion and is fi-finite, with No(5' > e) = 2/e for every e > 0. The notion 
of weak convergence is defined with respect to the finite measures No(-) = No(-,S' > e) (see e.g. [30]), 
and in particular n9\nt\ 7^^No(5 > t) = 2/{'jt). See [8, 39] for detailed surveys of super-processes 
and convergence towards them, and [11, 12, 36] for introductions to super-processes and continuous-state 
branching processes. 

In this paper, we study extensions of these results in the context of general spatial statistical mechanical 
models in sufficiently high dimensions that converge (or are conjectured to converge) to super-Brownian 
motion (SBM) in the sense of convergence of r-point functions. Convergence of r-point functions means 
that the (rescaled) joint moments of particle numbers and locations converge (to those of SBM). The use 
of r-point functions has a long history and tradition in statistical physics. The main result of this paper is 
that convergence of r-point functions, subject to two conditions that are valid in all our examples, implies 
that the classical results by Kolmogorov, Yaglom and (to some extent) Feller hold as well. As such, our 
result confirms that convergence of r-point functions is a relevant and important notion (see also [30] ) . 

Let us introduce the general setting that we investigate. Let P denote the probability measure de- 
scribing the law of our model. All our models have a notion of intrinsic distance, in which x — ^ y means 
that the shortest path between x and y has length n. Let Z+ = {0, 1, 2, ... } and R+ = [0, oo). Then for 
X £ and ft £ ZljT^ (or n G M^jT^ for models where time is continuous), we let 

tM(f) = p(0 ^ x,Vi = 1, . . . , r - 1) (1.1) 

denote the r-point function in the model. Further, for k = (ki, . . . , kr^i) G ([— tt, vr]'^)'"^"'^, we let 

t^;.\k)= ^ e^'^-%\x) (1.2) 

denote its Fourier transform, and 

e„ = P(3xGZ'^: O^x) (1.3) 

the survival probability. Let An = {x: x}, Nn = #{x: x}, and Sn = {Nn > 0} = {An 7^ 0}, 
so that On = P(5'„). When the underlying model is defined in discrete time, we define nt to be the vector 
{[nti\,. . . , [ntr\). 

In this paper, we investigate the asymptotics of the survival probability, assuming the asymptotic 
behavior of the r-point functions. These results apply to (a) lattice trees; (b) oriented percolation; and 
(c) the contact process, all above their (model-dependent) upper critical dimension, where the general 
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philosophy in statistical physics suggests that these models behave like branching random walk. In 
particular, when the allowed connections are sufficiently spread out, e.g. where all vertices within distance 
L ^ 1 of a vertex are considered to be neighbors of that vertex, the following condition holds as a theorem 
for each of these models, above their respective critical dimensions: 



Condition 1.1 (Convergence of the r-point functions), (a) There exist constants A^V > Q all depending 

,{r-l) 



on L such that for each r > 2 and t G 

1 



where the quantities (0) are the joint moments of the total mass at times ti, . 

measure of SBM. In particular, M*:~'>(0) = t'-'^2~^''-'^Xr - 1)!. 

tlr— 1 

(h) There exist constants A^V^v > all depending on L such that for each r > 2, t E 
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where the quantities Mj^ ^\k) are the Fourier transforms of the moment measures of the canonical mea- 
sure of SBM. 

Condition 1.1(a) is the weaker of the above conditions, and can be rephrased as 

r-1 

nE[ J] {Nt^n/n)] ^ A{VA^Y-^Mjr-'\0), (1.6) 

i=l 

where M^~^\0) are the limits of the joint moments of population sizes of critical branching processes 
with variance one offspring distributions. Note that the convergence in (1.6) makes no assumption on the 
spatial locations of the particles involved, however the evolution of Nn is affected by spatial interaction 
present in our models. Condition 1.1(b), which contains (a), can be rephrased as 
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r-1 




^Eno 











(1.7) 



where (j)k {x) = e for kj G M and x G Z , and where 



(1.8) 



Thus, Condition 1.1(b) states that certain moment measures of the rescaled processes under the measure 
/i„ converge to those of the canonical measure of SBM. Condition 1.1(b) is the condition that is typically 
proved in the literature. 

Before stating our main result, we start by formulating two further conditions. Let 



C(0) = {(x,n): ^ x} 



:i.9) 



denote the oriented cluster of G U^, i.e., all vertices x G Z'^ to which is connected, and we let |C(0)| 
denote its size. In continuous-time models, we instead take 



l-QO 

|C(0)| = / #{x: x}dt. 
Jo 

We make two central assumptions on our high-dimensional models: 



(1.10) 
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Condition 1.2 (Cluster tail bound). There exists a constant Cc such that 

P(|C(0)| > A;) < Cc/Vk. (1.11) 

Condition 1.3 (Self-repellent survival property). Let Tm be the a-field generated by the vertices at 
distance at most m from 0, i.e. by {{x,n): x,n < m}. Then there exists a constant Cg such that, 
with Nm equal to the number of x with x, almost surely for every stopping time M < n, 

F{Am n I Tm) < CeNM9n^M- (1-12) 

The cluster tail condition follows from the literature for all models under consideration. The self- 
repellent survival property in (1.12) turns out to be easy to check, and we shall do this below. The first 
of our main results is the following theorem: 

Theorem 1.4. When Conditions 1.1(a), 1.2 and 1.3 hold, as n ^ oo, 

nen^2/iAV), (1.13) 

and for each t > 0, 

f^n{Xi"\l) > 0) ^ No(Xt(l) > 0) = 2/t. (1.14) 

Consequently, conditionally on Nm > 0, the finite- dimensional distributions of {Nsn)s>t converge to those 
of Feller's branching diffusion started from an exponential random variable with mean A^Vt/2. 

For oriented percolation, the result reproves a result from [21, 22] (but without the error estimates) 
in a relatively simple way. See also [33, 34, 41] for related results on survival probabilities. Our set-up is 
rather general, so that in the future, it might be applicable to percolation and lattice animals as well. 

Theorem 1.4 is particularly important, since the combination of the convergence of the r-point func- 
tions (as formulated in Condition 1.1(b)) and Theorem 1.4 imply (see [30]) that {/x„}„>i converge in the 
sense of finite- dimensional distributions to Nq- This is the second of our main results. 

Theorem 1.5. When Conditions 1.1(b), 1.2 and 1.3 hold, the finite- dimensional distributions of the 
process (X^'"')(>o under fin converge to those of {Xt)t>o under the measure Nq. 

We now present some examples. All of the examples involve a function D\T!^^ [0,1], with 
SxgZ'* -^(^) — 1 t'^^'t obeys the properties of Assumption D in [26, Section 1.2] (whose precise form 
is not important for the present paper), together with [27, Equation (1.2)]. This assumption involves a 
parameter L G N, which serves to spread out the connections and which will be taken to be large. 



Spread-out oriented percolation above 4-1-1 dimensions. The spread-out oriented bond perco- 
lation model is defined as follows. Consider the graph with vertices Z'^ x and with directed bonds 
((x,n), (y,n -|- 1)), for n G Z+ and x,y G TL'^ . Let p G [0, H-Dll"^], where || • denotes the supremum 
norm, so that pD{x) < 1 for all x G Z"^. We associate to each directed bond ((x,n), {y,n -|- 1)) an inde- 
pendent random variable taking the value 1 with probability pD{y — x) and the value with probability 
1 — pD{y — x). We say that a bond is occupied when the corresponding random variable is 1 and vacant 
when it is 0. The joint probability distribution of the bond variables will be denoted by Pp, and the 
corresponding expectation by Ep. 

We say that {x,n) is connected to {y,m), and write (x,n) — > {y,m), if there is an oriented path 
from (x, n) to (y, m) consisting of occupied bonds. Note that this is only possible when m > n. By 
convention, (x, n) is connected to itself. We write {x, n) — > miim>n and there is a y G Z"^ such that 
{x,n) — > {y,m). The event {(0,0) — > oo} is the event that {(0,0) — > n} occurs for all n. There is 
a critical threshold pc > such that the event {(0,0) — > oo} has probability zero for p < pc and has 
positive probability for p > pc. The survival probability at time n is defined by 

^„(p) =Pp((0,0) ^n), (1.15) 
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and we let 6n = 0n{Pc)- General results of [3, 14] imply that lim„_!.oo Gn = 0. 

Then, for P = Fp^, Condition 1.1 is proved in [27]. Condition 1.2 holds by [2, 27, 37, 38], while 
Condition 1.3 follows from a union bound (i.e. F{Ux^Am{^ ~^ n}\!FM) < SxsAm ~^ ^I-^m)) and the 
strong Markov property. 

Spread-out contact process above 4 + 1 dimensions. We define the spread-out contact process as 
follows. Let Cn C be the set of infected individuals at time n G M_|_, and let Co = {0}. An infected site 
X recovers in a small time interval [n, n + e] with probability e + o(e) independently of n, where o(e) is a 
function that satisfies lim^^Qo{e)/e = 0. In other words, x G Cn recovers at rate 1. A healthy site x gets 
infected, depending on the status of its neighbors, at rate A J2yec„ ~ v)' where A > is the infection 
rate. We denote by P"** the associated probability measure. 

By an extension of the results in [3, 14] to the spread-out contact process, there exists a unique critical 
value Ac G (0, oo) such that 

0(A) = limP^(C„/0)(=°' (1.16) 

[ > 0, if A > Ac, 

and we define 

^n = ^n(Ac)=P^^(C„^0). (1.17) 

Condition 1.1 is proved in [23, 24]. Condition 1.2 holds by [23, 24, 1, 40], while Condition 1.3 again 
follows from a union bound and the strong Markov property. 

Spread-out lattice trees above 8 dimensions. A lattice tree is a finite connected set of lattice bonds 
(and their associated end vertices) containing no cycles. For fixed z > 0, every such tree T 3 with 
bond set B is assigned a weight Wz{T) = z'^' J|^^ y^^^ D{y — x), and we define Pz{x) = XItso x ^^ziT). 
The radius of convergence Zc of X^^g^d Pz{x) is finite. Let W{-) = WzX') aiid P = Pzd^)- We define a 
probability measure on the (countable) set of lattice trees containing the origin by P(r) = Given 
a lattice tree T 9 0, we define An{T) = {ai, . . . ,aN„} to be the (ordered) set of vertices in T of tree 
distance n G Z+ from the origin under some arbitrary but fixed ordering of U^. 

Condition 1.1 is the main result in [29]. Condition 1.2 follows from the detailed asymptotics for 
P(|T| = n) ~ cn^^/^ proved in [9, 10]. We next check Condition 1.3, for which it is enough to show that 
the result holds a.s. for every deterministic time m < n. Letting denote the tree up to tree distance 
m from the root, we have that F{Am — > n \ Tm = Tm) is equal to 

v-v '\/[/fT~'\ / > / J.i \^^/^{Ri avoid each other and Tm} {UjRj survives at least until n—m}i 

where X^/jg^ ^ swm over lattice trees R containing a G Z'^, and we recall that = {oi, • • • , otv^}- 

The final indicator function is bounded above by Ylj ^{SR.>n-m}-: where St is the survival time of T. 
By taking the sum over j outside and dropping the restriction that Rj avoids other Ri and t^, this is 
bounded above by 

E E w{R,)^sn.^>n-^^ 



(1.18) 



Y:T:T^=rW{T) 



i-i^j avoid each other and r„ 

<E E W{Rj)^{SR^>n-m} = NmpOn-m, 
j=l Rj^aj 
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where we have used the fact that the interaction term makes the graph Tm ^i^j Ri a lattice tree T with 
Tm = Tm, and weight W{T) = W {Tm)Y\i^jW {Ri) , so the numerator in brackets is no more than the 
denominator. This verifies Condition 1.3. 

Our main results can be restated in terms of the above models as follows: 

Theorem 1.6. Let L ^ 1, and let d > 4 for spread-out oriented percolation and the spread-out contact 
process, and d > 8 for spread-out lattice trees. Then, with A,V,v > all depending on L such that for 
each te M^"^^ and k e M^^-^) 



the asymptotics 

nen^2/{AV) and > 0) ^ No(Xt(l) > 0) = 2/t, as n ^ oo, (1.20) 

hold. As a consequence, the finite- dimensional distributions of the process {X'f^)t>o under /i„ converge 
to those of {Xt)t>o under the measure Nq. 

We close this section with two possible extensions to our results. 



Long-range models. In all our models, we assume that D has finite spatial variance, so that SBM 
can arise as the scaling limit. In the literature, long-range models have attracted considerable attention. 
See [5, 6, 7] for results on long-range oriented percolation, [18] for long-range self-avoiding walk, and 
[19] for percolation, self-avoiding walk and the Ising model. In long-range models, the random walk step 
distribution D has infinite variance. The simplest example arises when 



= . . e (1,21) 



where a G (0,2), and |x| denotes the Euclidean norm of x G Z'^. The results in [5, 6, 7] suggest that the 
upper critical dimension of oriented percolation equals 2a, while [19] indicates that it is 3a for percolation, 
and 2a for self-avoiding walk and the Ising model. 

We believe that Condition 1.1(a) holds for these models above their respective upper critical dimen- 
sions. Once this is proved, Theorem 1.4 then implies convergence of the survival probability in each 
case. However, random walk with step distribution D converges to a-stable motion rather than Brown- 
ian motion, a fact that is proved to hold for self-avoiding walk above 2a dimensions in [18]. Therefore, 
Condition 1.1(b) does not hold, and should be replaced with convergence towards the canonical measure 
of super-stable motion. 

By considering branching random walks, where the population size process is independent of the 
random walk step-distribution, it is easy to see that the law of the total mass process under the canonical 
measure of super-stable motion is the same as under Nq. Thus by [30, Theorem 2.6], in the long range 
setting, convergence of the r-point functions and the survival probability still implies convergence in the 
sense of finite-dimensional distributions. Therefore to prove a version of Theorem 1.5 in the long-range 
setting, it is sufficient to prove the convergence of the r-point functions in Condition 1.1(b). 

Spread-out percolation above 6 dimensions. Let p G [0, ||-D||oo^] be a parameter. We declare a 
bond {u,v} to be occupied with probability pD(v — u) and vacant with probability 1 — pD{v — u). The 
occupation status of all bonds are independent random variables. The law of the configuration of occupied 
bonds (at the critical percolation threshold) is denoted by Pp^ with corresponding expectation denoted 
by Ep^. Given a configuration we say that x is connected to y, and write x y, if there is a path of 
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occupied bonds from x to y, and the path with minimal number of bonds connecting x and y has precisely 
n edges. For percolation, Condition 1.1 is not known. The bound On < C/n is proved in [33] (in fact, 
we use Condition 1.3 together with an adaptation of the argument in [33] to prove that On < C/n in 
our general setting). Condition 1.2 follows from [15] together with [2], see also [16, 17]. As a result, for 
percolation, our results hold as soon as Condition 1.1 is proved. 

The above discussion suggests the following research program to identify the right constants in arm- 
probabilities in high-dimensional percolation, both in the intrinsic as well as in the Euclidean or extrinsic 
distance: (1) prove the convergence of the r-point functions in Condition 1.1(b) (from which the right 
constant in the survival probability or intrinsic one-arm probability would follow, improving upon the 
results in [33]); (2) prove tightness for convergence towards SBM; (3) identify the right constant for the 
extrinsic one-arm probability, improving upon the result in [34]. For the last step, an important ingredient 
showing that it is unlikely that a short path exists to the boundary of a Euclidean ball is proved in [25, 
Theorem 1.5]. 

The remainder of this paper is organized as follows. In Section 2, we prove an upper bound on On that 
is of the correct order, but with the wrong constant. In Section 3, we use weak-convergence arguments 
to identify the correct constant, and prove the consequences of convergence of the survival probability. 

2 Weak upper bound on the survival probability 

The following theorem gives a weak upper bound on the survival probability. 

Theorem 2.1. When Conditions 1.2 and 1.3 hold, there exists a constant c+ such that 

On < c+/n. (2.1) 

Proof. We follow [33], where a similar bound was proved for the intrinsic one-arm in percolation. We 
split O^n into two parts, 

Oin = ^{Nm >en\/mG [n, 3n], — > 4n) + F(3m G [n, 3n] : < en,0 — > An). (2.2) 

We can bound the first probability using (1.11), since |C(0)| > 2en^ if > en for ah m G [n,3n]. 
Therefore, 

C 

F{Nm >en'ime [n,3n],0 — > 4n) < P(|C(0)| > 2en^) < — ^. (2.3) 

nv2e 

In the second probability in (2.2), we let J > n be the first m G [n, 3n] such that < < en, and we 
condition on J"j = a{{Am)m<j)- Then, by (1.12), 

4n I Tj) < NjCeOn < enCgOn- (2.4) 

As a result, 

F{3m G [n, 3n] : Nn, < en,0 ^ 4n) = E[l|„< j<3„}P(^j 4n | Tj)] < eCguOl, (2.5) 
where we use the fact that n < J implies that — > n. Thus, we end up with the inequality 

Oin < + eCenOl. (2.6) 
nv2e 



Take e = C2 ^^'^ and take C2 > 1 so large that 



2-^2Cccj^ + Cecj^ <C2/A. (2.7) 

Then, it is easy to prove by induction that O^k < C24~^ for every k>l. By monotonicity of n 1— )• this 
immediately implies that On < (4c2)/n. This completes the proof of Theorem 2.1. ■ 
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3 Identifying the constant: Proof of Theorem 1.4 

In this section, we make use of general weak convergence arguments to prove that nOn — >■ 2/{AV). We 
rely crucially on a result that is essentially a special case of [30, Proposition 2.3], which requires the intro- 
duction of some more notation. Let Mi?(]R'^) (resp. Mi(]R'^)) denote the space of finite (resp. probability) 
measures on equipped with the topology of weak convergence. Let T^q denote the set of discontinuities 
of a function G, and D{E) denote the space of cadlag £'- valued functions with the Skorohod topology. 
When we say that ;U is a measure on (a topological space) E, this means that it is a measure with respect 
to the Borel u-algebra on E. 

Lemma 3.1. Suppose that Condition 1.1(a) holds. Then for every s,t,rj > 0, and every bounded Borel 
measurable H:R^R such that No(Xj(l) G Vh) = 0, 



f-i'TL 



^^(„)(^^^^^/7(X<"^(l))] ^ Eno [l{xAi)>v}H{Xt{M , asn^oc. (3.1) 



Proof. We follow the proof of [30, Proposition 2.3]. For convenience, we drop the superscripts 

By Condition 1.1(a), {/in}n>i is a sequence of finite measures on D{Mp{M.'^)) such that for every r > 1 

and t G [0, oo)'', (1.7) holds when cpkj = 1 for each j. 

Fix s,t,ry > 0. Let Ys = Xs{l), and define P„ = Pn,s,t G Mi{R'^) and P = Ps^t e Mi{m?) by 

. ^^Ays^{{Ys,Yt)&A}\ . _ IENo[:^.l{(y,,y06A}] 

where these measures are well defined since 

E^jy,] ^ENo[n] G (0,oo). 

On each of these spaces let {W,Z) be the canonical random vector, i.e. {W, Z){uji,uj2) = {uji,uj2)- Then, 
for every mi,m2 > 0, 

W rymi+lym2l jg rymi+lym2l 

Ep„ [^-iZ"^] = ^" V rv/ ^ ° V [V1 ' = > (3-2) 

^Mn [^s\ [ys\ 

i.e. the moments of (Wi Z) under P„ converge to those under P. 

Furthermore (see e.g. [30, Lemma 4.1]) there exists 5 > such that. 



E, 



< oo, (3.3) 



i.e. the moment generating function of {W, Z) under P is finite in a neighborhood of (0, 0). It then follows 
(see e.g. [4, Theorems 30.1 and 30.2, and Problems 30.5 and 30.6]) that P„ converges weakly to P, and 
therefore for G: ^ R bounded and such that P{{W, Z) G Vq) = 0, 

EpJG(iy,Z)] ^ Ep[G(VF,Z)]. 

In other words, for each bounded G: _^ M such that ^^{{Ys, Yt) G Vq) = 0, 

E^„ [YsG{Ys,Yt)] ^ Eno [YsG{Ys, Yt)] . 

Let H be as in the statement of the lemma, and define 



GH{x,y) 




\i X > r] 
otherwise. 
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Then Gh is bounded, and T>c^ = {(x,y): y G Vh or x = r/}, whence No((Xs,Xt) G T^Gh) — 0- The 
claim follows since ysGi^(ys5^t) = ^{Ys>-n}H{Yt). ■ 

Proof of Theorem 1.4- By Theorem 2.1, we have that n0„. is bounded. In order to investigate the limit of 
nOn, we split, for each fixed e > 0, 

nOn = nP(iV„ > en) + nP(0 < iV„ < en). (3.4) 

The first term is equal to {AV)^^ fin{X\"^ > ce), with c = (VA^)^^. From Lemma 3.1 with s = 1, rj = ce 
and with the continuous function H = I (and Condition 1.1(a)), we have that the first term on the right 
converges to (AF)-^No(Xi(l) > ce), and this converges to (^y)-iNo(Xi(l) > 0) = 2/{AV) as e 0. 
Since nP(0 < Nn < en) > 0, this immediately proves that 

lim inf n(9„ > 2/(^1^). (3.5) 

n— ^oo 

In order to identify the limit, we proceed as in [20]. Let 6 € (0, 1) and let {n^} = {nk{6)} be any 
subsequence of N such that rikOn^ — )• limsup„n^„ = 6, and (1 — 5)nk0(i_s)n^ — >• bs for some > 2/AV. 
This can be achieved by first taking a subsequence {ni/} for which miOmi b, and then taking a further 
subsequence {m/^.} such that (1 — 5)mi^6^i_^^^^ — )• hs- The required sequence is then n^ = mi^. 

Similarly to (3.4), for 5, e,e' G (0, 1) we write 

rikOn,, = nk¥{N(i_s)n^> £nk,Nn^> e'uk) 

+nfcP(iV(i_5)„^ > enfc,0 < iV„, < e'uk) + nfcP(0 < A^(i_5)„, < enfc,7V„^ > 0) 

= ^fc,<5,£,e' + Bk,5,e,£' + Dk,S,e- (3.6) 

Since the above is true for each 5,e,e', it follows that also 

b = lim sup rikOn^ < lim sup lim sup Ak^s^e^e' + lim sup lim sup Bk^s,£,£' + li™ sup lim sup Dk^s,ej (3-7) 

A;— >oo 5,e,e'4,0 fc— >oo 5, £,£'4-0 fc~>oo 5,£].0 k—^oo 

where the limits are taken in the order k — )• oo, e' | 0, e | 0, 5 | 0. 
The term Ak^s,e,e' can be rewritten as 

-l^^,^^{xt'J{l) > ce,xi"'=)(l) > ce') ^ -l^No(Xi_5(l) > ce,Xi(l) > ce'), as ^ oo, 

by Lemma 3.1. Letting e' J, and then e | this converges to 

-^No(Xi_5(l) > 0,Xi(l) > 0) = -l^No(Xi(l) > 0) = 2/AV, 

which, in particular, does not depend on 5. 

Further, using Condition 1.3, the term Dk^s^e satisfies 

Ce 

Dk,5,e = ^feIE[l{o<Ar(i_4)„^<£nfc}P(^^nfc > 0|J'(i_5)„j] < C0enk9snknkd(l-5)nk < _ 

uniformly in k, since nOn is bounded above uniformly in k. Letting e 4- 0, this converges to 0. 

We are left to investigate Bk^s,e,e'y ^or which we define, for each m, the measure Qm = IP(' I > 0). 
Then, we can rewrite 

^fc,<5,£,£' = nk9(i-5)rH:Q(l-5)rH:i^(l-5)nk > ^rik,0 < A^n^ < e'^fe)- 
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Thus, since nk9(^i_s)nk bounded above by < 2C for S < ^ (where C is independent of 6), proving 
that Umsup^^^ g/j^o l™sup;j_^oo Bk^s,£,e' = is equivalent to proving that 



hmsuphmsupQ(i_5)„j^(A^(i_5)„^ > enA:,0 < iV„* < e'"-fc) = 0. 

<5,e,e'4,0 k—^oo 

To prove (3.8), we note that, for any integers £1,^2 > such that ii + £2 > 1, 



1 





E 




-5)nk 




1 




\l-&)nk 




1 



I2 



(3.8) 



(3.9) 



n. 



where we use that N^i_^^^^ > when Nn^^ > 0, and where denotes a vector with precisely ii coordinates 
equal to (1 — 6)nk and £2 coordinates equal to n^. By Condition 1.1(a), 



n 



No 



Xl_5(l)^iXi(l 



(3.10) 



AV{1 - 5) 



No 



VA^Xi^s{l))''(yA^Xi{l)y'\Xi^s{l) > 



where the last equality follows from the fact that No(^i-<5(l) > 0) = 2/(1 — 6). Therefore, also using 
that (1 - 6)nk9(i-s)nf, ^ ^5, 



Er 



(l-i)nfe 



£2 



AVbs 



-E- 



No 



(vA^Xi.sil))'' [vA^Xiir 



Xi_5(l)>0 . 

(3.11) 



We recognize the above joint moments as the joint moments of {X, Y) with distribution (1— a5)5(o,o) +0^52^5, 
where (^(0,0) is the point measure on the vector (0, 0) and I's is the law of {A'^VXi^s{l), A'^VXi{l)) under 
No(-|Xi-5(l) > 0), and with as = 2/{AVbs). For any t > 1 - 5, 



,{Xt{l) = 0|Xi_5(l) > 0) = 1 - (1 - 5)/t, 



(3.12) 



so that 

us{Xi{l) = 0) = 1 - {1 - 6) = 6. (3.13) 

Let {Xn,Yn) be a two-dimensional distribution. Again by [4, Theorems 30.1 and 30.2, and Problems 
30.5 and 30.6], convergence of the joint moments of {Xn,Yn) to those of {X,Y) implies convergence in 
distribution when the moment generating function of both X and Y are finite in a neighborhood of 0. 
Under the conditional law No(-|-^i-5(l) > 0), the distribution of A'^VXi_s{l) is exponential with mean 
(1 - 6)A'^V/2 (see e.g., [20, Theorem 1.4]), and by (3.13), AVXi(1) is with probability S and an 
exponential with mean A'^V/2 with probability 1 — S. As a result, the distribution of both limits X and Y 
are mixtures of point masses at with probabilities 1 — as and l — as + asS and exponentials with positive 
means Ax and Ay. Therefore, their moment generating functions are finite in a neighborhood of zero, so 
that {N(^i_s)nk/nk, Nn^/rik) converges in distribution to (X, Y) having distribution (1 — a<5)(5(o,o) + oisfs- 
Thus, as fc — )• 00, 

Q(i_5)„,(iV(i_5)„, > enfc,iV„^ < e'rik) ^ asu^iA^V Xi_s{l) > e,A^VXi{l) < e'). 
When e' i 0, 

i/5(AVXi_5(1) > e, AVXi(1) < e') ^ MXisil) > ec,Xi{l) = 0) < ^^5(^1(1) = 0) = 5, (3.14) 
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where we use (3.13). Letting 5 | 0, we obtain (3.8). We conclude that Hm sup„_^oo = b< 2/{AV), 
which, together with (3.5), shows that liioin^oonOn = 2/{AV), as required. 

The fact that, conditionally on Nnt > 0, the finite-dimensional distributions of {Nsn/n)s>t converge to 
those of Feller's branching diffusion started from an exponential random variable with mean A'^Vt/2 can 
be obtained as follows. Fix t = sq < si < ■ ■ ■ < Sr < oo, and let < £ = J2]=o where £j, j = 0, . . . ,r 
are non-negative integers. Set s = (sq, . . . , Sr)- As in (3.9) we have 



n 

j=0 



n 



1 



■ n 



(0), 



tn 



where we now know that tnOtn — ^ 2/(Ay), and as before 

2 



n 



(0) 



AVt 



E 



No 



j=0 



Xt(l) > 



Thus the joint moments converge as in (3.11), i.e. 



n 

,j=0 



n 



II (V A' Xs^il)) ' Xt{l)>0 

j=0 



(3.15) 



Finally, the fact that {A'^VXs{l))s>t is Feller's branching diffusion follows from [13]. Again by [4, Theo- 
rems 30.1 and 30.2, and Problems 30.5 and 30.6], and the above bound on the moment generating function 
oi Xs{l), this completes the proof. ■ 
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